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The spin Seebeck effect is a spin- motive force generated by a temperature gradient in a ferromagnet 
that can be detected via normal metal contacts through the inverse spin Hall effect [K. Uchida et 
al., Nature 455, 778-781 (2008)]. We explain this effect by spin pumping at the contact that is 
proportional to the spin-mixing conductance of the interface, the inverse of a temperature-dependent 
magnetic coherence volume, and the difference between the magnon temperature in the ferromagnet 
and the electron temperature in the normal metal [D. J. Sanders and D. Walton, Phys. Rev. B 15, 
1489 (1977)]. 
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II. THERMALLY DRIVEN SPIN PUMPING 
CURRENT ACROSS f|N INTERFACE 



I. INTRODUCTION 

The emerging field called spin caloritronics addresses 
charge and heat flow in spin-polarized materials, struc- 
tures, and devices. Most thermoelectric phenomena can 
depend on spin, as discussed by many authors in Ref. [l|. 
A recent and not yet fully explained experiment^ is the 
spin anologue of the Seebeck effect — the spin Seebeck 
effect, in which a temperature gradient over a ferromag- 
net gives rise to an inverse spin Hall voltage signal in an 
attached Pt electrode. 

The Seebeck effect refers to the electrical cur- 
rent/voltage that is induced when a temperature bias 
is applied across a conductor. By connecting two con- 
ductors with different Seebeck coefficients electrically at 
one end at a certain temperature, a voltage can be be 
measured between the other two ends when kept at a 
different temperature. The spin counterpart of such a 
thermocouple is the spin current /accumulation that is 
induced by a temperature difference applied across a fer- 
romagnet, interpreting the two spin channels as the two 
"conductors". In Uchida et aVs experiment^ a temper- 
ature bias is applied over a strip of a ferromagnetic film. 
A thermally induced spin signal is measured by the volt- 
age induced by the inverse spin Hall effect (ISHE)^ in 
Pt contacts on top of the film in transverse direction (see 
Fig. [1} . This Hall voltage is found to be approximately a 
linear function (possibly a hyperbolic function with long 
decay length) of the position in longitudinal direction 
over a length of several millimeters. This result has been 
puzzling, since spin-dependent length scales are usually 
much smaller. The original explanation for this experi- 
ment has been based on the thermally induced spin ac- 
cumulation in terms of the spin thermocouple analogue 
mentioned above. However, the spin flip scattering shot- 
circuits the spin channels, and at which spin channels 
are shot-circuited, the signal should vanish on the scale 
of the spin-flip diffusion length^ 

In this paper, we propose an alternative mechanism in 
terms of spin pumping caused by the difference between 
the magnon temperature in the ferromagnetic film and 
the electron temperature (assumed equal to the phonon 
temperature) in the Pt contact. Such a temperature dif- 
ference can be generated by a temperature bias applied 
over the ferromagnetic filmj£ 

This paper is organized as follows: Section [TT] de- 
scribes how a DC spin current is pumped through a 
ferromagnet(F) [normal metal(N) interface by a difference 
between the magnon temperature in F and electron tem- 
perature in N. In section Mil we calculate the magnon 
temperature profile in F under a temperature bias. In 
Section HVl we compute the thermally driven spin current 
as a function of the position of the normal metal contact. 



In this Section we derive expressions for the spin cur- 
rent flowing through an F|N interface with a temperature 
difference as shown in Fig.[2J starting with the macrospin 
approximation in Subsection II. A and considering finite 
magnon dispersion in Subsection II. B. 

Since the relaxation times in the spin, phonon, and 
electron subsystems are much shorter than the spin- 
lattice relaxation time^^ the reservoirs become ther- 
malized internally before they equilibrate with each 
other. Therefore, we may assume that the phonon (p), 
conduction electron (e), and magnon (m) subsystems 
can be described by their local temperatures: yP' e ' m 
in F, and Tfj e in N4 We furthermore assume that the 
electron-phonon interaction is strong enough such that 
locally T£ = T e F =T F and = T%=T N . However, the 
magnon temperature may deviate: T™^Tp. This is il- 
lustrated below by the extreme case of the macrospin 
model, in which there is only one constant magnetic tem- 
perature, whereas the electron and phonon temperatures 
linearly interpolate between the reservoir temperatures 
Tl and Tr. The difference between magnon and elec- 
tron/phonon temperature therefore changes sign in the 
center of the sample. When considering ferromagnetic 
insulators, the conduction electron subsystem in F be- 
comes irrelevant. 



A. F N contact 

First, let us consider a structure such as shown in 
Fig. [21 in which the magnetization is a single domain and 
can be regarded as a macrospin M = M s Vm, where m 
is the unit vector parallel to the magnetization. We will 
derive in Subsection II. B the criteria for the macrospin 
regime. We assume uniaxial anisotropy along z, M s is 
the saturation magnetization and V is the total F vol- 
ume. The macrospin assumption will be relaxed below, 
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FIG. 1. (Color online) A ferromagnet (F) with thickness d 
and length L and magnetization M pointing in z-direction 
connects two normal metal (N) contacts at temperature T L / R 
at the left and right ends. A Pt strip of dimension Ixwxh 
on top of F converts an injected spin current (I s = I sp + 1/;) 
into an electrical current I c or Hall voltage Vh by the inverse 
spin Hall effect. 
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but serves to illustrate the basic physics. 

At finite temperature the magnetization order param- 
eter in F is thermally activated, i.e. m ^ 0. When we 
assume N to be an ideal reservoir, a spin current noise 
l sp is emitted into N due to spin pumping according to: 9 



At thermal equilibrium, T N = T™, and (I s ) — 0. At 
non-equilibrium situation, the spin current component 
polarized along m with prefactor in Eq. (JT]) averages 
to zero, thus does not cause observable effects on the DC 
properties in the present model. The x and y component 
of (I,) also vanish and: 



I«„(t) 



-in 



[g r m(t)xm(t) + gi m(i)] , 



(1) 



where g r and gi are the real and imaginary part of the 
spin- mixing conductance <7 m i x = 9r~^~^9i 

of the F|N inter- 
face. The thermally activated magnetization dynamics 
is determined by the magnon temperature T™, while the 
lattice and electron temperatures are Tp = Tp = Tp. 
The term proportional to g r in l sp has the same form 
as the magnetic damping phenomenology of the Landau- 
Lifshitz Gilbert equation (introduced below in Eq. ([5])). 
The energy loss due to the spin current represented by 
g r therefore increases the Gilbert damping constant. Ac- 
cording to the Fluctuation-Dissipation Theorem (FDT), 
the noise component of this spin-pumping induced cur- 
rent (from F to N) is accompanied by a fluctuating spin 
current Ifi from the normal metal bath (from N to F). 
The latter is caused by the thermal noise in N and its 
effect on F can be described by a random magnetic field 
h' acting on the magnetization^ 



7 



7 m(t)xh'(t). 



(2) 



In the classical limit (at high temperatures kBT^>fkoo, 
where ujq is the ferromagnetic resonance frequency) h' (t) 
satisfies the time correlation 

(7^(t)7*i(°)> = 2 -^^ L S ij S(t)=a' 2 5 ij S(t). (3) 

where (• • •) denotes the ensemble average, i,j = x,y, 
and a' = (^h/AiTM s V)g r is the magnetization damping 
contribution caused by the spin pumping. The correlator 
is proportional to the temperature T/y. 

The spin current flowing through the interface is given 
by the sum I s = I sp + I/; (see Fig. [5]). Here we are 
interested in the DC component: 



(Is) 



M S V 



[a' (mxm) + 7 (mxh')j , 



(4) 




FIG. 2. (Color online) F|N interface with 1) spin pumping 
current I sp driven by the thermal activation of the magneti- 
zation in F at temperature Tp, and 2) fluctuating spin current 
If 1 driven by the thermal activation of the electron spins in 
N at temperature Tjv- 
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[a'{m x rh y - m y m x ) - j(m x h' y - m y h' x }] 



(5) 

We therefore have to evaluate the correlators: 
(mi(0)rhj(0)) and (mi (0)^(0)). 

The motion of m is governed by the Landau-Lifshitz- 
Gilbert (LLG) equation: 



m = —7 mx (H c gz + h) + a mxrh, 



(0) 



where H g and a are the effective magnetic field and total 
magnetic damping, respectively, h accounts for the ran- 
dom fields associated with all sources of magnetic damp- 
ing, viz. thermal random field ho from the lattice asso- 
ciated with the bulk damping cto, random field h' from 
the N contact associated with enhanced damping a', and 
possibly other random fields caused by, e. g. additional 
contacts. Random fields from unrelated noise sources 
are statistically independent. The correlators of h are 
therefore additive and determined by the total magnetic 
damping a = a + a' + ■ • • : 

IrvyknT" 1 

<7/w(*hfci(o)> = ^ y f ¥(^MM- ( 7 ) 

The magnon temperature Tp is affected by the tem- 
peratures of and couplings to all subsystems: qT™ = 

a T F + a'T N H . 

We consider near-equilibrium situations, thus we may 
linearize the LLG equation. To first order in mx<l 
(with m z ~l) 



+ ariiy = -uj m y + -fh y , 
- am x = +uj m x - jh Xl 



(8a) 
(8b) 



where ujq = "fH c g is the ferromagnetic resonance (FR) 
frequency. With the Fourier transform into frequency 
space g = Jge wt dt and the inverse transform g — 
fge~ luJt duj/2ir, Eq. © reads m,-(a)) = J2jXij(u)"fhj(u), 
with i,j = x, y and the transverse dynamic magnetic 
susceptibility 



xM = 



wq — ictoj —iu) 



(ujq - iaijj) 2 - up- \ iu Wo - iauJ 
in terms of which, utilizing Eqs. (|31[7]), 



(9) 



K(f>,(o)) = — 



cr 2 /"XtjM-XjiH _^ t du 



(m i (t)h' :j (0)) = — jxiMe-^ — . 



(10b) 
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Eq. (|10ap gives the mean square deviation of m 
in the equal time limit t — > 0: (mj(O)mj(O)) = 
SitfkBTp/woMsV. The time derivative of Eq. (jTUal) is 



2ir 



By inserting Eqs. (llObt [TTT) with t into Eq. ©: 
M,V /a' 



(11) 
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2 /A f r , , / \ -I dUJ 

~ a - a J \-Xxy(U) -Xyx(u)\ 



1 + a 



B (rpm rp \ jft-grks (rprn rp \ 



2ttM s V 



N 



(12) 



where L' s is an interfacial spin Seebeck coefficient. 
In Eq. (fT2"|) . we used Jim. [x^ (w)] du/2w — (since 
Eq. (jlObl) is real and Inr^ changes sign when uj —± —cj), 
and / \Xxy(u) - Xyx(u)] (dw/2%) = 1/(1 + a 2 ) (see Ap- 
pendix [K§ . From Eq. (fT2")l we conclude that the DC spin 
pumping current is proportional to the temperature dif- 
ference between the magnon and electron/lattice temper- 
atures and polarized along the average magnetization. 

When the magnon temperature is higher (lower) than 
the lattice temperature, the DC spin pumping current 
flows from F into N leading to a loss (gain) of angular 
momentum that is accompanied by the heat current: 



Q, 



(Iz)Hctf — K' m A(Tp l — Tjv), (13) 



where A is the contact area and K' m = 
u)Q[j,BkB(gr/A)/irM s V is the interface magnetic heat 
conductance with Bohr magneton /ib- 

In addition to the DC component of the spin pumping 
current, there is also an AC contribution to the frequency 
power spectrum of spin current and spin Hall signal. 11 
A measurement of the noise power spectrum should be 
interesting for insulating ferromagnets for which the large 
imaginary part of the mixing conductance can be much 
larger than the real part (see Appendix |B|) . 



B. Magnons 

In extended ferromagnetic layers the macrospin model 
breaks down and we have to consider magnon excita- 
tions at all wave vectors. The space-time magnetization 
autocorrelation function can be derived from the LLG 
equation (see Appendix [Cjl : 



(rhi(0,0)mi(0,0)) 



M s V a 



(14) 



where we introduced the temperature-dependent mag- 
netic coherence volume 



V a 



3Z & 



4nD 



(15) 



with D the spin stiffness, and Z the Zeta function. Phys- 
ically, this coherence volume, or its cube root, the coher- 
ence length, reflects the finite stiffness of the magnetic 
systems that limits the range at which a given perturba- 
tion is felt. When this length is small, a random field has 
a larger effect on a smaller magnetic volume. 

The results obtained in Subsection II. A can be car- 
ried over simply by replacing V — > V a in Eqs. (J3J- 
[7J. The corresponding spin Seebeck coefficient is L' s — 
-fhk B {g r /A)/2TTM s V a , and the heat conductance is K' m = 
uJol J -BkB(gr/A)/T:M s V a . Here we assumed that the 
magnon temperature does not change appreciably in the 
volume Va^V. 



III. MAGNON-PHONON TEMPERATURE 
DIFFERENCE PROFILE 

Sanders and Walton (SW) e discussed a scenario in 
which a magnon-phonon temperature difference arises 
when a constant heat flow (or a temperature gradient) 
is applied over an F insulator with special attention to 
the ferrimagnet yttrium iron garnet (YIG). Its antiferro- 
magnetic component is small and will be disregarded in 
the following. This material is especially interesting be- 
cause of its small Gilbert damping, which translates into 
a long length scale of persistence of a non-equilibrium 
state between the magnetic and lattice systems. SW as- 
sumes that, at the boundaries heat can only penetrate 
through the phonon subsystem, whereas magnons can 
not communicate with the non-magnetic heat baths (i.e. 
K' m = in our notation). Inside F bulk magnons in- 
teract with phonons and become gradually thermalized 
with increasing distance from the interface. The different 
boundary conditions for phonons and magnons lead to 
different phonon and magnon temperature profiles within 
F. However, according to Eq. (TT31 . the magnons in F 
are not completely insulated as assumed by SW when 
the hear reservoirs are normal metals. In this section, 
we follow SW and calculate the phonon-magnon temper- 
ature difference in an F insulator film induced by the 
temperature bias but consider also the magnon thermal 
conductivity K' m of the interfaces. In case of a metallic 
ferromagnet the conduction electron system provides an 
additional parallel channel for the heat current. 

As argued above, the boundary conditions employed 
by SW need to be modified when the heat baths con- 
nected to the ferromagnet are metals. In that case, the 
magnons are not fully confined to the ferromagnet, since 
a spin current can be pumped into or extracted out of 
the normal metal. In Fig. [TJ the two ends of F are at dif- 
ferent temperatures: Tl and Tr of the N contacts drive 
the heat flow. Let us ignore the Pt contact on top of 
the F for now. When both phonons and magnons are in 
contact with the reservoirs, by energy conservation (sim- 
ilar to Ref. @) the integrated heat Q mp {Qmp) flowing 
from the phonon to the magnon subsystem in the range 
of —L/2 < z' < z (z < z' < L/2) has the following form: 
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imp 



imp 



(z) 



CpCrn 1 



' mp 



[T p {z')-T m {z')]dz' (16a) 



-K n 



dT m 

1 dz 



K' m [T m (-L/2)-T L ] 



-K P ^-K' p [T p ^L/2)-T L ], 



CpC m 1 



' mp 



1 dz 
dT r , 



[T p (z') - T m (z')]dz' 

Z 

K' m [T m (L/2)-T R ] 



-K p ^-K' p [T p (L/2)-T R ], 



(16b) 
(16c) 

(16d) 

(16e) 
(16f) 



where C p>m (Ct = C p + C m ) are the specific heats, K pm 
(Kt = K p + K m ) are the bulk thermal conductivities for 
the phonon and magnon subsystems, K p m (K' T = K' p + 
K' m ) are the respective boundary thermal conductivities. 
T mp is the magnon-phonon thermalization (or spin-lattice 
relaxation) timeii The boundary conditions for T m and 
T p are set by letting z = ±L/2 in Eq. (fl6|). i.e. 



K 



dT, 



m.p 



m,p 



dz 



z=±4 



TK' mtP [T mtP (±L/2)-T R/L ], (17) 



The solution to Eq. ([TBI and Eq. (|17p yields the magnon- 
phonon temperature difference AT mp (z) = T m (z) — 
Tp(z): 



AT mp (z) 



K T {K' p K m - K' m K p ) sinh f AT 



-K m K p {K' T L + 2K T ) cosh ^ + {2K' m Kj + 2K' p K 2 m + K' p K' m LK T ) sinh ^ 



sinh j 
sinh^ 



AT 



(18) 



with AT = Tl — T R and 



X 2 = 



C v + C n 



K p K n 



CpCrn 



K p + K n 



I m .p r 



Cf, 



(19) 



where the approximation applies when C p ^$>C m and 
Kp^K m . Eq. ()18|) shows that the deviation of the 
magnon temperature from the lattice (phonon) tempera- 
ture is proportional to the applied temperature bias and 
decays to zero far from the boundaries with characteristic 
(magnon diffusion) length A. 

We use the diffusion limited magnon thermal conduc- 
tivity and specific heat calculated by a simple kinetic 
theory (assuming hujo^ksT)^ 



C r n. — 



15C(|) 
32 



k%T 3 

7T 3 T> 3 



and K„ 



mi) 



16 



k 7 T 5 t 
7T 3 T> H? 



with T m the magnon scattering time. Using these expres- 
sions in the approximate form of Eq. (|19p we obtain: 



A = 



14C(7/2) Dk B T_ 
3C(5/2) h 2 



(20) 



In Appendix [Dl we estimate T mp ~l / (2auJo) for ferromag- 
netic insulators, assuming that magnetic damping a is 
caused by magnon-phonon scattering. It is difficult to 
estimate or measure r m , and the values quoted in the lit- 
eratures ranges from 10 _9 s to 10 _7 s, depending on both 
material and temperature^*^ At present we cannot pre- 
dict how A varies with temperature: Eq. (|20[) seems to 
increase with temperature, but the relaxation times likely 
decrease with T. 



When KL 



and K p 



— > oo, i.e. 



the boundary is 



thermally insulated for magnons and has zero thermal 
resistivity for phonons, the prefactor r/ in Eq. (|18p reduces 
to SWs result^ 



V - 



A 



K T 

Xpcothi + 2K, 



m 



l 



(21) 



AT, 



nip 



where the approximation is valid when K p ~^>K. 
is obviously maximal in this limit. 

The discussion in this section also applies to ferromag- 
netic metals when the electron-phonon relaxation is much 
faster than the magnon-phonon relaxation. The electron 
and phonon subsystem are then thermalized with each 
other and can be treated as one subsystem. In this case, 
we may replace K p — > K pe = K p + K e and 



n n 1 

^p^m 1 

Cp ~\~ C m T m p 



CpC m 



1 



CpC„ 



1 



Cp ~h Cm 



C e ~\~ Cm T me 



(22) 



IV. HALL VOLTAGE 

In Uchida et al.'s experiment (Ref. 0), a Pt contact is 
attached on top of a Py film (see Fig. [JJ to detect the 
spin current signals by the inverse spin Hall effect. A 
spin current polarized in the z-direction that flows into 
the contact in the y-direction is converted into an electric 
Hall current I c and thus a Hall voltage Vh — I C R in the 
x-direction. 
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YIG 


Py 


Unit 


7 


1.76xlO n 


1. 76x10" 


1/T-s 


4ttM s 


a 1.4xl0 5 


/ 8.0xl0 5 


A/m 


D 


b 1.55xl0- 38 


/ 7.6xl0~ 39 


J-m 2 


a 


a 5xl0~ 5 


s 0.01 


— 


LU 


no 


s 20 


GHz 


Tmp 


c,d w -6 


d 10- 7 (Ni) 


s 


Tm 


c,e^Q — 9^7 


no- 9 


s 


Qr 
A 




no 18 


1/m 2 


K 1/3 


5.4 


3.8 


nm 


V 


0.4 - 0.5 


0.27 


— 


A (th) 


4.7 - 47 


0.3 


mm 


A (exp) 


a 6.7 


n.o 


mm 


t (th) 


0.38 - 3.8 


130 


/iV/K 


£ (exp) 


a 0.16 


fc 0.25 


txV/K 



TABLE I. Parameters for YIG and Pyjat T = 300K if not 
specified, LT « 10 K). a Ref. E3, 6 Ref. LLS c Ref. GJ, d Ref. 0, 
e Ref. Q, 'Ref. M {Dis derived from A ex = 13 pJ/m), 9 Ref. 
M, h Ref. M, *Ref. M 3 Ref. 



Quantity 


Values 


Reference 


6h 




0.0037 


4 


P 




0.91 pfi-m 


,2 


Ixwxh 


4 mm 


x 0.1 mm x 15 nm 


2 



TABLE II. Parameters for Pt contact. 



For the setup shown in Fig. [T] we assume that the Pt 
contact is small enough to not disturb the system, which 
is valid when the heat flowing into Pt is much less than 
the heat exchange between the magnons and phonons or 

K' m AT mp (z)(lxw)^ C, f m 1 AT mp (z){lxwxd), 

which is well satisfied when d~^>l nm for both YIG and 
Py at low temperatures. 

From Eq. (fP8"]l . we see that at the position below the 
contact (z — z c ) the magnon temperature deviates from 
the lattice (and electron) temperatures by T™ — T F = 
AT mp (z c ) . If we assume that the contact is at thermal 
equilibrium with the lattice (and electrons) in the F film 
underneath, i.e. T c e = TJP = T F (z c ), then a temperature 
difference between the magnons in F and the electrons 
in Pt exists: Tf{z c ) - T c e = AT mp (z c ). Therefore, by 
Eq. (H2J), a DC spin-pumping current (I z ) from F to Pt 
is driven by this temperature difference, which gives rise 
to a DC Hall current in the Pt contact (see Fig. [I): 



h A 



-zxy, 



(23) 



where 6h is the Hall angle. In Eq. (|2"3"|) we disregarded the 
spin diffusion backflow and finite thickness corrections to 



the Hall effect, which is reasonable when thickness h is 
comparable to the spin diffusion length i s( j. 

According to Eq. (|12[) , the thermally driven spin pump- 
ing current can flow from F to Pt: 

(I z ) = L' s [T m (z c ) - T c ] = L' s [T m (z c ) - T p {z c )]. (24) 

Making use of Eq. (fl"8|) , the electric voltage over the two 
transverse ends of the Pt contact separated by distance 
I is V H (z c ) = pljc(z c ): 



Vi 



H 



V^Hplgrieks sinh- 
irM s V a A sinh 



-AT=£ 



sinh ■ 



2A 



sinh 



-AT, (25) 



2A 



Using the numbers in Table U and Table [XT] and 
a value of A~7 mm that reflects the low magnetiza- 
tion damping in YIG, we have ^~0.47 from Eq. (|2ip . 
and £~0.38~3.8AiV/K from Eq. {25]) for YIG for 
g r /yl~10 15 ~ 16 /m 2 , which is consistent with with the ex- 
periments by Uchida et al. (Ref. [TH). For Py, we estimate 
£~130 [iV/K, which is almost three orders of magnitude 
larger than the experimental value of 0.25 /iV/K given 
in Ref. U 

The crucial length scale A is determined by two ther- 
malization times: the magnon-magnon thermalization 
time r m and the magnon-phonon thermalization time T mp 
as shown in Eq. (|2Tj|) . Knowledge of these two times is 
essential in estimating A. The quoted values of r m mp 
are rough order of magnitude estimates at low temper- 
atures. Yet, in order to completely pin down the value 
of A for this material, a more accurate determination of 
T m ,mp and K m , C m as a function of temperature by both 
theory and experiments is required. As seen in Eq. (|25p . 
the spatial variation of the Hall voltage over the F strip 
is determined by the magnon-phonon temperature differ- 
ence profile as calculated in the previous section. From 
Eq. (f2"D)) and the parameters for YIG in Table Q] (where 
the values of r m and T mp are very uncertain), we estimate 
A~4.7~47 mm, which is again consistent with the exper- 
imental value of A~6.7 mm. For Py, we estimate A~0.3 
mm (using the r mp value for Ni instead of Py), which is 
about one order of magnitude smaller than its measured 
value. 



V. DISCUSSION & CONCLUSION 

A magnon-phonon temperature difference drives a DC 
spin current, which can be detected by the inverse spin 
Hall effect or other techniques. This effect can, vice 
versa, be used to measure the magnon-phonon tempera- 
ture difference. Because the spatial dependence of this ef- 
fect relies on various thermalization times between quasi- 
particles, which are usually difficult to measure or calcu- 
late, this effect also accesses these thermalization times. 
On a basic level, we predict that the spin Seebeck effect 
is caused by the non-equilibrium between magnon and 
phonon systems that is excited by a temperature basis 
over a ferromagnet. Since the inverse spin Hall effect 
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only provides indirect evidence, it would be interesting 
to measure the presumed magnon-phonon temperature 
difference profile by other means. Such measurements 
would also give insight into relaxation times that are dif- 
ficult to obtain otherwise. 

In principle, the theory holds for both ferromagnetic 
insulators and metals. However, as shown above, the 
agreement between the theory and experiments is rea- 
sonably good for ferromagnetic insulator YIG, but the 
theory fails for ferromagnetic metal Py, which underes- 
timates the length scale A and overestimates the mag- 
nitude £. This might be because of several reasons: (i) 
we have completely ignored the short-circuiting effect of 
the metallic Py, to which the inverse spin Hall current 
leaks, (ii) the lack of reliable information about relax- 
ation times T mp>m for Py ccould cause the difference in 
A, (iii) the complication due to the existence of conduc- 
tion electrons in ferromagnetic metals. 

In conclusion, we propose a mechanism for the spin 
Seebeck effect based on the combination of: (i) the in- 
verse spin Hall effect, which converts the spin current 
into an electrical voltage, (ii) thermally activated spin 
pumping at the F|N interface driven by the phonon- 
magnon temperature difference, and (iii) the phonon- 
magnon temperature difference profile induced by the 
temperature bias applied over a ferromagnetic film. Ef- 
fect (ii) also introduces an additional magnon contributed 
thermal conductivity of F N interfaces. The theory holds 
for both ferromagnetic metals and insulators. The agree- 
ment between experiments and theory is satisfactory for 
insulating ferromagnet. The magnitude and the spatial 
length scale for YIG is predicted to be in microvolt and 
millimeter range. The lack of agreement for both the 
length scale and the magnitude of the spin Seebeck effect 
for permalloy remains to be explained. 
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vanishes by Jordan's Lemma): 



Xe 



2tt 



-i Im 



1 



-Ci-a v 



©(<)• (Al) 



The minus sign comes from the counter-clockwise con- 
tour, and 0(i) is the Heaviside step function, which van- 
ishes when t < and is unity otherwise. This integral 
is discontinuous at t = 0, therefore the value at t = is 
given by the average of the values at t — 0^: 



-^o^ = 1 1 
X 2tt 2 1 + a 2 




(A2) 



Appendix B: Spin-mixing conductance for F|N 
interface 

In this Appendix, we use a simple parabolic band 
model to estimate the spin pumping at an F|N inter- 
face, where F can be a conductor or an insulator. The 
Fermi energy in N is Ep, and the bottom of the conduc- 
tion band for spin up and spin down electrons in F are 
at Ep + Uq and Ep + Uq + A, respectively, where Uq is 
the bottom of the majority band and A is the exchange 
splitting. The reflection coefficient for an electron spin a 
(t or 4-) from N at the Fermi energy reads 



r(*) = 



k -h k n 



(Bl) 



y/2moEp /h 2 — q 2 is the longitudi- 
in N (q is the transverse wave- 
vector), fcj- = y/2m,f(EF — Uo)/h 2 — q 2 and 



where k = 
nal wave-vector 



k± = 



^2 mi {E F - U 
vectors (or imag 
spins. m is the 
fective mass for 
reads 



A)/h 2 — q 2 are the longitudinal wave- 
inary decay constants) in F for both 
effective mass in N and m a is the ef- 
spin <7 in F. The mixing conductance 



9n 



A f 

— / (l - r t rj) d 2 q = g r + %g t . 



(B2) 



Appendix A: Integral 



We evaluate Eq. (IB2[) with m = m a having the free 
electron mass, E F = 2.5 eV, U = - 4 eV.— A 
plot of the mixing conductance is shown in Fig. [3] for 
A = 0.3,0.6,0.9 eV. Uq = corresponds to a ferromag- 
netic metal (with majority band matching the N elec- 
tronic structure), and Uq > Ep = 2.5 eV corresponds 
to a ferromagnetic insulator (with \r a \ = 1 for both spin 
types). 



Here we evaluate the frequency integral in Eq. (fT2"|) . 
To this end, we need to reintroduce the time dependence 
and then take the limit t — > 0. When t > 0, the in- 
tegral in Eq. (jlObp can be calculated using contour in- 
tegration (real axis + semi-circle in the lower plane so 
that e~ lujt — > 0, where the integral over the semi-circle 



Appendix C: Magnetization correlation for 
macroscopic samples 

The dynamics of m is governed by the LLG equation: 
m(r,t) = -7m(r,t)xH ef t + am(r, i)xrh(r, f), (CI) 
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/ - gi :A = 0.3eV - V 




] g r : A = 0.6 eV \ 




- gi : A = 0.6 eV 




g r :A = 0.9eV 




-gi:A-0.9eV 





1 2 3 

barrier height: U (eV) 



FIG. 3. (Color online) g r ,i vs. Uo in logscale for different A 
values. Curves starting with values close to unity give the 
real part g r , and curves starting from very low values display 
the imaginary part — gi. The vertical axis is in units of the 
Sharvin conductance. 



where H cff = H + (D/jh) V 2 m(r, t) + h(r, t) is the total 
effective magnetic field with: (i) the external field plus 
the uniaxial anisotropy field Ho = Hqz = (uio/-f)z, (ii) 
the exchange field (D/'yh)V 2 m due to spatial variation of 
magnetization, and (iii) the thermal random fields h(r, t). 
We define Fourier transforms: 



g(K u) = I dre^J die*" g(r, t) (C2) 
^ = vE e-* kr fpe^g(k,u) (C3) 

V * — *m } n,l J Z7T 
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where k = 2n(m/L,n/W,l/H) and L,W,H are the 
length, width, and height of the ferromagnetic film. Near 
thermal equilibrium, m±<^l (and so m z c±V), we may 
linearize the LLG equation. After Fourier transforma- 
tion, the linearized LLG equation becomes mj(k, uS) — 
HjXij(^^)lhj{\s.,u) with i,j = x,y, and 



X = 



1/(1 + a 2 ) ( ojk — iaoj —iuj 

(bJ — LJ^ ) (bJ — 0J^ ) I ILU bJfr ~ ictOJ 



(04) 



with fiwk = hujQ + -D|k| and = ±Wk/(l±ia). 

The spectrum of random motion in the linear response 
regime Xi = J^jXijfj with "current" x and random force 
/ (x and / are chosen such that xf is in the units of 
energy) is comprehensively studied in Ref. H|| In our 
problem, / ->• jh, x -> (M s /j)m, \ ~> ( M s/l)x, the 
autocorrelation of the magnetization then becomes (£ = 
ri -r 2 ,r = ti - t 2 ) 



(mi(CT)mj(0,0)) = 



2 1 k B T 



M S V 

iujrXij ~ 



Xji h\oj\/k B T 



2ir 



I-0J 



^h\uj\/k B T 



(C5) 



from which (rhi(C T)rrij(0, 0)) can be caluclated by tak- 
ing derivative over r. 

The limit C = an d r = can be obtained for three- 
dimensional systems by replacing the sum over k by an 
integral V -1 ^ / d 3 k/(2n) 3 in Eq. (IC5l) : 



K(o,ok-(o,o)) = 8 tj , 



(^(0,0)771^(0,0)) = 



3Z^k B T /fc R r 



(k B T\ 



2M S 



(C6) 

-a -l\ 
1 -a) ' 
(C7) 



where ^ is the Zeta function. 



Appendix D: Relation between a and r mp 

Here we derive a relationship between the magnetic 
damping constant a and the magnon-phonon thermal- 
ization time r mp . When the magnon and phonon tem- 
peratures are T m and T pi the magnon-phonon relaxation 
time is phenomenologically defined by£ 



d T —T 
— (T - T ) - --^ — 



(Dl) 



If the phonon system is attached to a huge reservoir (sub- 
strate) such that its temperature T p is fixed (dT p /dt = 0), 
Eq. (jDip becomes 



T — T 



' mp 



(D2) 



In metals we should consider three subsystems (magnon, 
phonon, electron) leading to 



dt 



-^k st (T s -T t ), 



(D3) 



with s, t = m, p, e for magnon, phonon, and electron. 



The coupling strength k s 



l C t /(C s + C t ) with the s- 



t relaxation time T s t and the specific heat C s . When the 
magnon specific heat is much less than that of phonons 
and electrons (C m <cC p , C e ) 



dT m 
dt 



T — T T — T 



(D4) 



' mp 



We may parameterize the magnon temperature by the 
thermal suppression of the average magnetization: 



fleffM a V(l - (m z )) = k B T m , 



(D5) 



where H c q points in the z-direction. The LLG equation 
provides us with the information about m z : 



m = — 7mx(F£ ff + h) + a mxm, 



(D6) 
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with a p / e the magnetic damping caused by scattering 
with phonons / electrons and a — a p + a e . Therefore 



where the random thermal field h = h p + h e from the lat- 
tice and electrons are determined by the phonon/electron 
temperature: 

K/.(t)7Aj /e (0)> = 2l ^' M k f^ (D7) 
i I 

k B dT m 1 2~fk B (aT m - a p T p - a e T e ) 

- (m z ) = {m v -fh x - m x jh y ) + a{m x m y - m v m x ) = - ■ ^ -^-^ , (D8) 



H cS M s V dt v y 1 * * 1 yl x y y 1 1 + a 2 M S V 

r 



or wo-10 GHz (Ref. [Tj, thus T mp = l/2aw -10" 6 s, which 



For the ferromagnetic insulator YIG: a~6.7xl0 and 

mp 



dT m 2cjq 

~T7~ = — ; — ^{aT m ~a p T p -a e T e ), (D9) 
at 1 + or 



grees with r mp ~10 6 s in Ref. 



ae 



with uiq = ■yHcff. Comparing Eq. (ID9[) with Eq. (|D4[) . we 

find The estimate above relies on the macrospin approxima- 

2 2 2 1 1 tion. Considering magnons from all k adds a prefactor 

2au Q ~ aUJ ° = apUJ °, + aeUJ °, = + . (D10) 01 OTder 01 unity 

1 + a z 1 + a z 1 + a 1 T mp r me 
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